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SUMMARY 

The longitudinal stability characteristics of elastic swept 
wings of high aspect ratio experiencing bending and torsional 
deformations are calculated for '"supersonic speed by the appli- 
cation of linearized lifting-surface theory. A parabolic wing 
deflection curve is assumed and the analysis is simplified by 
a number of structural approximations . The method is thereby 
limited in application to wings of high aspect ratio for which the 
root effects are small. Expressions for the lift, pitching-moment, 
and span load distribution characteristics are derived in terms of 
the elastic properties of the wing; namely, the design stress, the 
modulus of elasticity, the shearing modulus, and the maximum design 
load factor. The analysis applies to wings with leading edges swept 
behind the Mach lines. In all cases, however, the trailing edge is 
sonic or supersonic. Application of the method of analysis to wings 
with leading edges swept ahead of the Mach lines is discussed. 

The results of numerical calculations for a wing of aspect 
ratio 3.2 and 60° sweepback axe presented for a Mach number of 
lAl^ and for incompressible flow. The effects of wing elasticity 
on the lift— curve slope, moment— curve slope, and neutral— point 
position are shown. The results indicate that the primary variable 
involved in aeroelastic phenomena is the dynamic pressure and that 
the influence of the flight Mach number is small for wings swept 
behind the Mach lines. 


INTRODUCTION 

In reference 1, R. T. Jones has shown that supersonic flight 
may be attained with a reasonable degree of efficiency through 
the use of swept wings of high aspect ratio. The use of sweep- 
back, however, involves many problems of stability and control, 
not the least of which are associated with the aerodynamic effects 
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of the elastic deformation of the airplane structure. In particular, 
the longitudinal stability cf the aircraft may he affected to a large 
degree since the tending and torsional deformations of the wing may 
shift the center of pressure cf the lift forward an appreciable 
distance . 

These aeroelastic phenomena occur under those flight conditions 
where the magnitude and/or the spanwise variation of the elastic 
deformation of the wing varies with angle of attack. Aeroelastic 
effects may therefore occur either in accelerated flight at constant 
dynamic pressure or, under certain conditions, in steady level flight 
with varying dynamic pressure. In the latter case, if the loading 
due to twist or camber is different than the loading due to change 
of angle of attack, the trim change due to elastic deformation of the 
wing in steady level flight varies" with the dynamic pressure and 
influences the stability of the airplane as indicated by the position 
of the control stick as a function of airspeed . 1 

In solving aeroelastic problems, since the interrelation of the 
structural and aerodynamic characteristics of the wing results in 
mathematical complexity, it is usually necessary to compromise to 
some extent either the structural or the aerodynamic aspects of the 
problem to obtain a solution. In the present analysis, the structural 
characteristics of the wing are compromised to the extent that the 
form of the deflection curve is assumed. This assumption permits the 
application of supersonic lifting— surface theory to the determination 
of the load distribution, the lift, and the pitchi n g-moment character- 
istics of elastic wings. Additional analysis is necessary to deter- 
mine whether it is better to use more rigorous aerodynamic theory in 
aeroelastic computations, as in the present report, or to use a more 
complete structural theory as in recent work by John W. Miles of 
U.C.L.A. or Franklin Diederich of the HACA. 


SIMBOLS 

Xi,yi Cartesian coordinates 

x,y tr ans formed Cartesian coordinates in terms of the semispan 

dimension, s 


1 This partic ular aeroelastic characteristic is not considered in the 
present report which is concerned primarily with accelerated flight. 
Further, the wing is considered to be weightless so that the 
ameliorat ing influence of the distributed mass of the wing is not 
accounted for in estimating the aeroelastic characteristics. 
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x,y coordinates of the apex of any superposed lift ing 
sector 


s distance in the y x direction from the root section to the 

intersection of the flexural axis and the tip Mach cone 

s* distance along the flexural axis from the root section to 

the intersection of the flexural axis and the tip Mach 
cone 

y* distance measured from the root section along the flexural 

axis 


Y 


S 

X 


c 


c 


AH 


spanwise distance in y direction from the root section to 
the center of load on the half wing 

wing area 

taper ratio, ratio of tip chord to root chord 
average chord 


mean aerodynamic chord 


f I c 2 dy \ 
\ / c dy ) 


local chord parallel to the plane of symmetry 

root chord parallel to the plane of symmetry in terms of 
the span dimension, s 

aspect ratio 


A angle of sweepback of the flexural axis 

6 slope of the flexural axis in a vertical plane passing 

through the flexural axis 

n maximum load factor 


M bending moment at any point on the flexural axis 




bending moment at the root section of the w ing beam 


T torsional, moment at any point on the flexural axis 
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?r 

E 

G 

I 

J 

£ 

°max 

*r 

a 

O&p 


torsional moment at tlie root section of the wing "beam 
modulus of elasticity for the wing beam material 
shearing modulus for the wing beam material 
moment of inertia of the wing team 
torsional stiffness constant 

distance between the flexural axis and the center of 
pressure of the sectional lift in terms of the local 
chord 

maxi mum design stress 

nw--sH imim thickness of the wing at the root section 

angle of attack of the root section of the wing 

incremental angle of attack at any spanwise station of the 
wing 


ctg angle of attack of the wing section at any spanwise station 

a angle of attack of the root section at which maxi mu m load 

n factor is developed 

3 ,Jyp — 1 where M is the free— stream Mach number 

m 3 t ime s the cotangent of the angle of sweeptack of the 

wing leading edge 

nij. p times the cotangent of the angle of sweeptack of the 

wing trailing edge 

t p times the cotangent of the angle of sweeptack of a ray 

from the apex of any superposed lifting sector 

E complete elliptic integral of the second kind with modulus 

(*/i=52) 

W airplane weight 


a 



ra|*J 
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<L 

AP 

q. 

1 


wing loading 


dynamic pressure f ipY 2 ^, where p is the mass density and 
V the velocity a? the free stream 


lifting pressure coefficient 
load per unit span 


c^ section lift coefficient 


L lift 

C L lift coefficient 

lift coefficient at maximum load factor 

mg section pitching moment of a wing section about the apex 

of' the wing 



C m pitching-moment coefficient about the apex of the wing in 

terms of the mean aerodynamic chord and the wing area 

C-j-^ the rate of change of lift coefficient with the angle of 
attack of the root section 


C-m^ the rate of change of pitching-moment coefficient with the 
angle of attack of the root section 

Cbv, the rate of change of pitching-moment coefficient with the 
I* lift coefficient 


ANALYSIS 

Wing With a Subsonic Leading Edge 

In the following analysis, for convenience, the aerodynamic 
loading due to bending and that due to torsion are first treated 
separately. Expressions for the combined effects of bending and 
torsion are derived later. 

Bending.— The aerodynamic twist 2 due to bending of a stream- 
wise section of an elastic swept wing under accelerated flight 
2 

The c hang e in camber of the airfoil sections due to the distortion 
of the wing surface Is, of course, ignored. 
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conditions is a function of the applied load and the elastic character- 
istics of the w ing team. In order to arrive at a solution for the 
aerodynamic properties of the wing without becoming involved in labo- 
rious graphical analysis, some simplifying approximations must be made 
regarding the elastic properties of the wing. 

In a strict sense, a swept wing of conventional structural design 
cannot be considered to have a flexural axis. For wings of high 
aspect ratio, however, it will be assumed that a flexural axis exists, 
since this assumption permits the use of simple beam theory and intro- 
duces only a small conservative error. 

For the purpose of analysis, the root section of the wing beam 
is assumed to be the extension of the wing beam on a plane perpendic- 
ular to the flexural axis and passing through the intersection of the 
flexural axis and the streamwise root section. (See fig. 1.) This 
simplification of the beam analysis is similar to that of reference 2. 
The length of the wing beam s' is taken as the distance along the 
flexural axis from the root to the intersection of the flexural axis 
and the tip Mach cone. The semispan s of the wing is taken as 
extending from the root section to the intersection of the flexural 
axis and the tip Mach cone in a direction perpendicular to the plane 
of symmetry. The portion of the wing lying within the tip Mach cone 
is ignored since, as shown in reference 3 S very little load is carried 
in this region and the analysis is thereby simplified. 

The coordinate system is selected as shown in figure 2. The 
origin of the coordinate system is placed at the apex of the wing, 
the positive branch of the Xi axis lying downstream. 


The mathematical treatment may be made Issb tedious by trans- 
forming and nondimens ionaliz ing the coord in ates so that in the follow- 
ing analysis 


7 = 


s 


X 


Zl 

s 


_ _ root chord 

c o s 


In general, at both subsonic and supersonic speeds, selection 
of the wing plan form for low drag leads to a combination of spanwise 
loading and spanwise distribution of the "bending resistance in the wing 
"beam such that the wing deflection curve is essentially parabolic. ( The 
ratio of M to I is constant across the span.) The deflection curve 
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deviates appreciably from a parabola only if the aeroelastic effects 
experienced by the wing are very large. Calculations show that, for 
the usual tapered -wing, the assumption of a parabolic deflection 
curve gives results comparable to a more rigorous structural trea tment 
for deviations from rigid wing values as large as, for instance, a 
30-percent loss in lift— curve slope. It seems improbable that a 
designer would be interested in wings with larger deviations from' 
rigid-^rtng characteristics . 


Since the deflection curve of the flexural axis is assumed to be 
parabolic, the slope of the flexural axis is 


e 


= JL yt 
El 


where y* is measured along the flexural axis. 


The incremental angle of attack of streamwise sections of the 
wing is related to the slope of the flexural axis as 


ocrp = — 0 sin A 

The slope of the flexural axis in nondimensional transformed coordi- 
nates may be written as 

M s 

9 y 

El p cos A 

The incremental angle of attack of any streamwise section of the 
elastic wing is then 

Ms, 

“T = ~ M P 7 t£m A 

and the total angle of attack of any streamwise section is 

M s 

os = a - — - y tan A (l) 

where a is the angle of attack of the root section of the wing. 
Equation (l) gives the magnitude and distribution of twist across 
the span of the wi n g if the magnitude of M/EI is known. 

The distribution of pressure over the elastic wing due to twist 
may be determined by applying known conical-flow solutions for super- 
sonic flow. In the linearized theory, the principle of superposition 
of various solutions may be used to satisfy the particular boundary 
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conditions of the problem. For the elastic ving, the flov field 
may be considered to consist of the superposition of tvo distinct 
flov fields: 

1. The flov about a flat rigid ving at an angle of attack 

equal to the angle of attack of the root section 

2. The flov about a tvisted ving for vhich the angle of 

attack at the root is zero 

The solution for the first flov field is given in references 
4 and 5) the second flov field can be obtained by determining the 
solution for a differential tvist da s at one station and inte- 
grating this solution across the span. 

The solution for the pressure distribution at any point, corre- 
sponding to a differential tvist, must meet the foil owing boundary 
conditions (fig. 3): 

1. Outboard of the station of tvist, the angle of attack must 
be constant end equal to the differential tvist. 

2. Inboard of the station of tvist, the angle of attack of the 
surface must be zero. 

3. Between the swept leading edge and the Mach cone, no lift- 
ing pressures may exist. 

The conical— flow solution corresponding to these boundary condi- 
tions is that for a special lifting sector given by Lagerstrom in 
reference 6 and is expressed in the notation of the present report as 

Ad _ 8a m 3//2 / 1+t 

q j5lt m+1 v mr-t (2) 

vhere t defines a ray from the apex of the sector. 

Figure 3 shows both a sketch of the boun da ry conditions to be 
met by this solution and a plot of the pressure distribution given 
by equation (2). 

The induced pressure resulting from tvist due to bending of the 
elastic ving may be found by integrating across the span of the ving. 
This integration corresponds to the superposition of an infinite number 
of the lift ing sectors along the span, each sector having an infinites- 
imal angle of attack da s . 
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The pressure due to twist is then given "by 


where 



8 m 3/2 Ms 
p 2 * (nH-l) ^ 


tan A 



dri 


t = 2=1 = m (y~T)) 

x— £ mx— T] 

The x and y coordinates of the apex of any superposed sector 
are | , *1 . 

The Integration must he carried out from the root section cf the 
wing t\=0 to the value of 11=110 corresponding to the laBt superposed 
sector, the Mach cone of which encompasses the point x,y under con- 
sideration. The value of tj 0 is found hy placing t equal to —1 and 
solving for tj. 


‘ ^ (l+y) 


The integration yields at any point . x,y the pressure due to 


twist 



16 m 5/a , .Ms 

— 5 ; To tan A — 

3P 2 * (mKL) 2 El 


(x+y) H** 

hi mx— ; 


(3) 


To this expression must he added the conjugate term due to the elastic 
deformation of the opposite wing panel. The conjugate term may he 
obtained hy substituting — y for y. Then 



16 m 5 / 2 , ■ . Ms 

tan A — 

3P 2 * (mH) 2 El 


(x+y) + (x-y) 

L J mx-y J mxfy J 


(4) 


It should he noted that the addition of the conjugate terms adds 
some very small lifting pressure in the region between the wing leading 
edge and the Mach cone where no lifting pressure may exist. These 
pressures may he canceled hy the superposition of constant lift sectors 
as noted in reference 7* Since these extraneous pressures, on the 
average, amount to about 3 percent of the average pressure coefficient 
over the adjacent wing surface and, since elimination of those pressures 
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would change the pressures over the surface only about one— half of 
1 percent. It seems that in view of the additional complication 
involved the cancellation of these pressures is unwarranted. 

The total lifting pressure for the elastic wing at an angle of 
attack is then obtained by adding to equation ( 4 ) the solution for 
the flat lifting wing. For the elastic wing, then 

( Z +y) 

V mx— y 


Ap kmPa l6 

BE J m 2 _(Z) 2 3 P 2 * (mhl) 2 


m 5 / 2 , * Ms 

tan A — 

EE 


+ <*-*> m - 


(5) 


Examination of this equation shows that the relationship between M/EE 
anil a must be established before the pressure distribution can be 
calculated. Since for wings with parabolic deflection curves the 
nwHimim stress occurs at the point of maximum thickness, usually 
the root, the maximum stress occurring at maxi mum load factor is 


q max ~ 



2 


ana since the bending moment at any point on the span is a linear 
function of the angle of attack. 


M _ 2cfTnax CL 
I dp CLq 

where Umax is the design stress at m a ximum load factor, d r is the 
thickness of the root section and an is the angle of attack at 
maximum load factor; an expression for oq is derived later. 

The equation for the pressure distribution may then be written as 


Ap _ 4 m g a 


(2+y) /S + (x-y) /S' 


(6) 


3P 2 Jt (m+-l) 2 


E < 3 -r °tn 


mx-y 


mx+y J 
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The load per unit span can he obtained from an integration of equation 
(6) with respect to 2 along any streamwise station (yi=constant) 


2 

4 





m 



d2 


The integration is carried out from the leading edge of the wing, 
x = ^ to the trailing edge 2 = and yields 

*11 ~m 1 


1 

q. 


4m 2 sa 

PS 


fx(y) 


-iA d_ ten A 283 £ 2 . f 2 (y) 

3P 2 * (mfl) E ^ On 


(7) 


The functions fi(y) and f 2 (y) are given in the appendix since they 
are somewhat unwieldy. 

The lift coefficient may he obtained by an integration of equar- 
tion (7), spanwise from root to tip. 


P C L 

The integration yields 3 



bc =2sf[i!s5f:F 32 m 5/2 
L S L BZ 1 ^B 2 * (m+l} 


The constants F x and F 2 are given by equations in the appendix. 

This equation may be used to determine the angle of attack at 
maximum load factor which is needed in the foregoing equations: 

®It may be noted that the ratio s 2 /S is essentially the same as 
one— fourth aspect ratio and that the parameter s/d r is directly 
related to s f /d r , a common structural criterion. 



12 


mCA TN No. 1811 


^ ~ 4m2Fi 


s 

2b2 


°Ln + 


32 

33 s * 


m 5 / 2 

(m+l ) 2 


Umax 

tan A 

E 


s 

d r 



( 9 ) 


The pitching-moment characteristics of the elastic wing may he deter- 
mined. hy an integration of the pressure distribution given by equation 
( 6 ) . 


For any spanwise station, the section pitching moment about the 
apex of the wing is 


^b=-b 2 


This integration yields 


' ”t (Si) xto 

V Q. >y^const. 


m 


^ = -s 2 

a 


5/2 


knPa , . 32 m , u max 

% ( 7 ) - tanA 


s a , . I 
U (j) 

dr c&n J 


pZ - - - 3 p*n Utt+l}" 1 E d r eta 

The functions £3 (y) and £* (y) are given in the appendix. 


( 10 ) 


The total pitching-moment coefficient about the apex of the wing 
in terms of the mean aerodynamic chord is found by integration across 
the span. 


pc 


m 


_ 2s 
Sc J 


2 s 


[ 1<m2a _ 32 


5/2 


or 


ftp — 2 s 3 

s <r 


L & 

pz 


m 

dy 
cl ^ 


tan A JE££ JL ±. jv 
E dr Qn J 


3p2jt (mH)2 

5/2 ■ 

32 m . . %ax s ^ 4 


(n) 


( 12 ) 


The constants E 3 and E 4 which are functions of the aspect ratio, 
taper, and sweepback are given in the appendix. 

Torsion .— The previous analysis has ignored the effects of wing 
twist due to torsion. The solutions obtained are, in reality, those 
for wings of infinite torsional stiffness. In general, since the 
flexural axis (or torsion center) is behind the center of pressure at 
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all spamd.se stations of the wing, the twist of the wing due to 
torsion will tend to compensate for the twist due to henfling. For 
wings having large angles of sweepback such as are necessary for 
efficient supersonic flight, the aerodynamic twist due to torsion 
has been calculated to he about 15 to 20 percent of the twist due to 
bending (for thin wings). In such cases, the effect of the torsional 
deformation on the spanwise loading may be neglected in calculating 
the torsional moment. Equation (7) may be utilized in the calcula- 
tion of the torsional moment in this instance. A complex simulta- 
neous solution is thereby avoided. 

An expression for the torsional moment at the root section of 
the wing beam (perpendicular to the elastic axis) may be obtained by 
assuming that the distance from the center of pressure to the flex- 
ural axis for any section of the wing is a constant percentage of the 
local chord. 

Then 

T 

= s ^ (cos A) / -c dy 
J o 1 

where c is the local streamwlse chord given by the equation 


c 


s c 0 


1 — (l— a) 

P. 


Where X denotes the taper ratio of the wing and t the distance 
from center of pressure of flexural axis in terms of the streamwise 
chord. The function describing the spanwise loading l/q is given 
by equation ( 7 ). 


as 


The equation for the torsional moment at the root may be written 



£ Co cos 



1 - (1— x) | 



b 2 i 


Co cos A 



Idy-i (1-*) 





Ab will be shown later, it Is convenient to derive the ratio of the 
torsional mome nt at the root to the bending mo m ent at the root. The 
bending moment at the root Is given as 



Ik 
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and 



P cos A 



£r 


c Q cos 2 A 





( 13 ) 


and corresponds to the spanwise center of pressure for the load on the 
h a lf wing. The value of T may he determined by a mechanical or 
analytical Integration of equation ( 7 ). 


When the assumption is made that the twist due to torsion varies 
linearly across the span (or that the ratio T/GJ is constant across 
the span), the incremental angle of attack of any section of the wing 
due to torsional deflection may he written as 

_ cos A Ts T „ Ts y 
^ 0 GJ 7 * GJ 0 


or 



and hy adding this expression to the angle of twiBt due to bending 
(equation (l) ) the total angle of twist of any section is 


“e " a ” Ef "p tanA 

-(SS) 

T r ' 

J 


(lAa) 

<x e - a. 2 = °»y * I [ 

nr E On P L 

tan A — 

(S) 

1 1 

(lkb) 


Combined bending and torsion .— Expressions for the aerodynamic 
properties of swept wings experiencing both bending and torsional 
deformation may be obtained from equations (6) to (12) if tan A is 
replaced by 


tan A — 
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The equation for the angle of attach at m a x i Tmrm load factor for 
combined “bending and torsion is then 


4m S F 1 

15^ 0L “ 

where 

Tr 


Mr 


+ 32 m5/g _s_ 
30 3 * (mil) 2 E d r 


tan A — 


Ji'G' Mp -J -J 


(15) 


= i c Q 


L f " (1 ^> 


COB 


2 A 


In applying the foregoing analysis to a specific ving, it is 
convenient to use the equations to obtain the ratio of or 

Cj^ for the elastic wing to the value for the rigid wing. Multi- 
plying this ratio by the value of or for the rigid wing 

as determined by the complete theory wherein the region within the 
Mach cone of the tip, and so forth, is considered, will then give 
more accurate parameters for the elastic ving. Then 


l 

qa 


4m 2 s _ / x _ 32 m 5/g g max s^ f g {y) 
£E 1 (mil) 2 E dp 


tan A — 


JpG^Mr 


(16) 




’‘'elastic 1 _ 8 E Jm. 

Gl «rigid 3 P«(mH) 2 


max s 


E dpCji 


tan A — 


L^XT-plEg 

JrG/Mr J F x 


(IT) 




’elastic _ 


>rigid 


1 _ 8 E ^/m g max s 
3 Pit(m+l) 2 E dpO^ 


tan A — 


IySATrlEt 

J r G/M r J E3 


(18) 


In using the preceding equations, it Is necessary to solve for a^. 
This, in turn, involves finding the ratio T^My which is determined 
by the parameter T (usually has a value of about 0. 4o) . 


A solution of the combined bending and tors_ional deformation 
effects can be obtained by assuming a value of T, solving for a^, 
and checking the value of T from a moment anri area integration 
of a plot of equation (l6) to see if a second approximation is required 
to determine more accurately. 
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The previous equations apply primarily to flat lifting wings or 
to twisted and cambered vings for which the loading due to twist and 
camber is the same essentially as the loading due to change in angle 
of attack. For wings with somewhat arbitrary camber and/or twist, 
these equations apply to all accelerated flight conditions. A solution 
for the aeroelastic characteristics in steady level flight for such 
wings must involve a consideration of the effects of the loading due 
to the known arbitrary twist. 


Wing With Supersonic Leading Edge 


The foregoing analysis has treated wings with the leading edge 
swept behind the Mach cone. The same method, however, may be applied 
to wings swept ahead of the Mach cone. In this case, however, the 
expression for the pressure field for the incremental twist at any 
spanwise station, corresponding to equation ( 2 ) , is given by refer— ; 
ence 8 as the real part of 


Afi = ifcS gk cos -1 fHpfc 

q pir^/mS-l | t-m | 

where <x,t, and m are as defined for equation (2) . 


(19) 


Expressions for the pressure distribution, lift, moment, and load 
distribution may be obtained in the same manner as for a wing with a 
subsonic leading edge although the integrations are more involved. 


DISCUSSION 

Supersonic Lifting-Surface Theory 

The results of the foregoing analysis are best illustrated by 
applying them to a specific wing. For this purpose, the wing shown 
in figure 4 was selected, having the geometric and structural material 
characteristics given in the table in the figure. The calculations 
were made for various values of the parameter nh and for two values 
of the maximum design stress. 4 ® 

Span load distributions for the wing are shown in figure 5 for 
a Mach number of 1.4-14, a value of n^ of 150 pounds per square foot, 

s 

4 Calculations show that the wing has sufficient depth to withstand the 
maximum loading assumed without failure. 
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a design stress of 30,000 pounds per square inch, and a dynamic 
pressure of 211 pounds per square foot -which corresponds to flight 
at 60,000 feet altitude. The load distribution curves of part (a) 
of figure 5 are for the same angle of attack of the root section and 
show that the elasticity of the wing results in an appreciable 
decrease in lift— curve slope. In this case, the reduction experienced 
by the elastic wing amounts to 15 percent of the value for the rigid 
wing of the same plan form. Part (b) of figure 5 shows the load 
distribution curves for constant total lift coefficient. These load 
distributions are of significance in illustrating how the change in 
span load distribution due to elasticity may be expected to shift 
the longitudinal center of pressure forward. The load distributions 
as derived by what is known as strip theory are discussed later. 


Comparison of Aeroelastic Effects at Supersonic Speed 
With Incompressible Plow Solutions 

In calculating lift and stability characteristics of elastic 
wings, it should be noted that errors resulting from assuming the 
extent of the wing beam as given in figure 1 and from ignoring the 
lift within the tip Mach cone may be minimized by using the analyti- 
cal expressions which give the ratio of lift-curve slope or the 
ratio of moment-curve slope for the elastic wing to that for the 
rigid wing. These ratios may be used with the rigorous values of 
and from reference 3 to obtain accurate values of 

and for the elastic wing. 

Such ratios have been computed for the wing shown in figure 4 
as functions of the dynamic pressure at a flight Mach number of 1.4l4. 
For comparison, the same ratios have been computed as functions of 
the dynamic pressure for incompressible flow by the theory of refer- 
ence 9* Figures 6 and 7 show the results of these calculations 
which were made for two value b of nj| of 150 and 300 pounds per 
square foot and two values of design stress, 30,000 and 45,000 pounds 
per square inch. Figure 8 shows the shift in neutral point 5 due to 
wing elasticity as calculated from the data of figures 6 and 7. 

The results indicate that the aeroelastic phenomena are a little 
more severe at supersonic speed, although the aeroelastic effects are 
found to be primarily a function of the dynamic pressure and not of 
Mach number . At a given dynamic pressure the differences In the 
aeroelastic effects aB computed by incompressible flow theory and by 

Neutral point is defined as the position of the center of gravity 
along the mean aerodynamic chord for neutral stability. 
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supersonic lifting— surface theory are found to be due largely to the 
fact that the center of pressure of the sectional lift is farther 
forward at subsonic speed, resulting in a difference in torsional 
deformation which compensates somewhat for the tending deformation. 
The comparison indicates that the dynamic pressure is the primary 
■variable involved in determining the aeroelastic characteristics, 
at least for wings swept behind the Mach lines. 

It should be noted that the variation of the aeroelastic charac- 
teristics with dy nami c pressure given in figures 6, and 8 for any 
value of and o vra-sr includes the effect of a small variation in 

wing beam moment of inertia which comes about from the manner in 
which the rrenrlmirm design Btress was brought into the analysis. The 
effect is not significant within the range of dynamic pressure for 
which the theory applies. 

Tn regard to the range of application of the equations, calcu- 
lations made using more rigorous structural theory with simple strip' 
theory show tha t the method of the present report may be expected 
to give accurate estimates of aeroelastic effects as great as, for 
instance, a 30 — percent loss in lift-curve slope. Within such limits 
it is expected that the estimate of the neutral point shift due to 
elasticity will be much more accurate than for analyses using 
elementary aerodynamic loading. 


Strip Theory 


The analytical evaluation of aeroelastic effects can be greatly 
simplified by the use of strip theory. This simplified method of 
determining the aerodynamic loading is based on the assumption that 
the loading at any spanwise station of the wing is a function only 
of the section angle of attack. For the present case, a modified 
strip theory suggests itself wherein only the incremental lift change 
due to elasticity is considered to be a function of the incremental 
local angle— of-attack change due to elasticity. The loading at any 
spanwise station is then given by the product of the ratio of the 
local angle of attack to the angle of attack of the root section and 
the expression for the rigid wing loading at any spanwise station. 
Then, for strip theory 


qa, 


a 


The load distributions so calculated are compared with supersonic 
lifting— surface theory on figure 5. It is evident that this form 
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of strip theory overestimates the effects of wing elasticity, hut 
the comparison indicates that the accuracy of strip theory in 
predicting lift-curve slope is satisfactory. The shift in center of 
pressure which strip theory gives, however, is much too conservative. 

It is suggested, however, that the modified form of strip theory 
may prove very useful in estimating the effects of wing elasticity on 
certain aerodynamic parameters, the damping in roll for instance. 
Further, the use of modified strip theory permits the structural 
characteristics of the airplane wing to he brought into the problem 
more conqpletely and enables the designers to estimate wing character- 
istics for all modes of deflection. 


Ames Aeronautical Laboratory, 

national Advisory Committee for Aeronautics, 
Moffett Field, Calif. 


APPENDIX 

MATHEMATICAL EERIVATION OF LOADING FUNCTIONS 
AND PLAN-FCEM CONSTANTS 

The functions f x (y), fa(y)^ % (y), and ft (y) and the constants 
F 2 , F^ y ^4 which appear in equations (7) to (18) of the text are 
given in this appendix. 


The functions fi(y) and f 2 (y) were developed from the 
following integral: 


l 

q = S 




i” (■ 


from which 


y 

m 


q 'y^const. 


dx 


y™t c o 
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■which yields 


fi(y) 




y2 + c D y + m2c 0 2 


f2(y) = / 


which yields 


y-®tc 0 


y-HUfcCo 


(w) dx + / -0 (x-y) 


mx+y 


f 2 ( ) = 2 m(y-Httbc 0 )-mty(^m+ 3 ) 


[y-Hntc 0 +mty) (my-Hnmfcc 0 -mtyJ 


" ly 2 (m+l) 2 l T . . -x.-i 2 m(y+mbC 0 )-«afcy(m-l) 

L 'J L COSl1 mfc(m + l)y J 


+ ) 7 (y+m-tfo^ty^ (my+nmi fc0o +m t yy _ 


[ 3y a ( m+ l) g ] [ cosh -i gm(y-fmfcCo)^y(in^l) 1 

L 8m®/ 2 J L mt(m+l)y -* 

f y 2 (^m +l) y 2 (l-m ) *1 [ SrKm+l ) 2 [ . _i 3 -m 

+ L km. 5 / z J L 8m 5 / 2 J L m+1 


The constants Fj. and F 2 are evaluated as follows : 


F3. = / fi(y) dy 
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which yields for mj. ^ 1 


n = sw5®- ! 


2m-^ma 


2a 


, m tf 2 ^ _! m^-ap _! m > \ 

+ — ( cos x cos ± 


2a 3 / 2 


mrn-^f 


n^y 


and for m-^ = 1 


mc f 


*1 = \4a 2 Co 2 +2pm 2 c 0 +p 2 (m 2 — l) 

2m 2 ( 1-m 2 ) 2 (1-^n 2 ) 


Co" 


2 (1-m 2 ) 3 / 2 


cos 1 


,-1 _ cos-1 m 


me. 


where 


a = m^ 2 ~m 2 


f = 




and 


*2 = P f 2 (y)dy 

^o 


which yields for mt / 1 


F * - (jssps) [ (rf2)3/2 - <^ 2 hA) 3/2 ' 
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£^1] 


mt 2 (mfl)^f 2 
8b ^ 


_, df— 2f3b _! d 

COS 1 — < COS A — r- 

mtf (bH- 1 ) mt (nH-1 ) 


} 


’ p 3 (mfl) 2 " 


L 8m 5 /^ J 



' f (m+1) 2 ' 

\ 

' f 2 (3d 2 +4mb ) 

r coo -i ,.-1 4 ..] 

8m 2 

[ 8b s/ 2 

mtf (mfl ) mt (m+1 ) 


3<lf+2pb 

“ ~ F“ 


Vmf 2 +pdf— p 2 ^ 


3m 1/2 df £ 

“4b 1 


M 


p 3 (5nri-l) 72 (1-m) 


12m5/2 



” (mf 2 + phf — p 2 g ) 3 / 2 — (mf 2 ) 3 / 2 

3g 




+ 


f (4mg-4ij) l / 2pg-faf 
8m 2 mt 2 g \ 4g 


//mf 2 +phf— P 2 g + 


m 1 / 2 ^ % 

4g 


+ 


f 2 mt 2 (1+m) 

8g 3 7 2 


2 r 


_ n bf-2pg _! h 

COS ; r COS — r 

mtf (m+l) mt (m+1) 


} 


f (mfl) 2 ] / f 2 (31i 2 +4mg| 

8m^ J l 8^ 


cos 


1 


bf— 2Pg 
mtf (mfl) 


cos 


1 


h 

mtCm+i) 
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where 


For mt = 1 


b = (mt-ai) (mt+1) 
d = 2m+mmt-mt 
e = 2m+5mmt+3 m t 
f = mtc 0 
g = (m+mt) (mt - l) 
h = 2nMrmt+mt 
j = 3mt+5™it-2m 


F 2 = f f[mco 2 +Pco( 3 b+- l)+2p 2 (m— 1) ] 3/a _ (mc 0 2 ) 3/2 \ 

\ 4m 2 J l 6(nt-l) 6(m— l) J 


£o _ c 0 (3m-l)(7m+3) 

2m l6m 2 (m-l) 


{ 


r 4g(mr-l)+c 0 (3a^-l) 
8(m-l) 


Vmc 0 2 +pc 0 ( 3m— 1) +2p 2 (m-1) 


c 0 (3m-l) 

. 8(m-l) J 


(*/mc 0 2) 


c 0 2 (m+l) 2 


Ll 6 3 / 2 J 


cos-i ^P(ni-l)+Co(3m-l) 
c 0 (m+l) 


— cos 


-i 3m-l 


m+1 J 


P 3 (m+l) 2 

8m 5 / 2 


} 

cosh- 1 M.3?a- 1 )+^ c 
p(m+l) 


c 0 2 (m+l) 2 
8m 

3c 0 (3m-l) 

_ l6m 1/ ' 2 (m-l) 2 _ 


• l|p(m-l)-3c 0 (3m-l ) 
l6mc 0 (m— l) 2 


cosh- 1 _ cosh - 

P(nH-l) 


*/mco 2 +Pco( 3®— 1) +2p 2 (m— 1) 


-i 3—™. 1 

m+1 J 



2k 
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c n (l9m g — lQm+3) 

L 32 v'Sa(l-m) 5/2 J 


cos 


1 


^fi(m-l)+c 0 (3nfr~l) 

c 0 (m+l) 


cos 


3m-l 

m+1 _ 


f2 [SmSco 2 — M3mc 0 (m+l)+33 2 (m+l) 2 1 1 

- 

i6m 3/,2 c 0 

\zkc 0 +P(m+l) 

15 mc 0 3 / 2 (m+l ) 3 J 

15 (m+1) 3 
» 


jr i j l 0 ^o 0 (»i) T /a } 

ll_3m(m+l)J |_lQm 2 c 0 (m+l) JJ LL J J 


[■ 0 3 ( 5 m+l) \/ 2 (l— m) 

8 (mc 0 2 .) 3/2 1 

12 m 5 / 2 J 

_ 15 m(m+l) J 


The functions f 3 (y) and f 4 (y) were developed from the integral: 


2§_ _ 

q 


y+mt c o 

m t 


= — s‘ 


/. ( 


^ /y=const. 


m 


from which 


i(7) = x 


y + mt c o 

“t 


xdx 


J m 2 — (y/x) 2 


m 


Then 


f 3 (7) = ^ m2 < 7 ' +mtCo) ' 2 " mt ^ 
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y 2 m(y+^tc 0 ) 

+ — n cosh. *■ 

2m 0 mj-y 


Also 


y-^gp 
mt 


y +m t c o 


^ = x u+t) y^ iaz+ x <**> / 


x-y 

mx+y 


m 


m 


which yields 


f4(y) = 


y(y+mtCo+mty) (my+mmtco-mty) 


8 (y +m fc c o) J 

2km mt 3 


2mm t y(y-Hn t c 0 ) ( 7m+5 ) +m fc 2 y 2 ( 3 m 2 + 22 m+ 15 ) 

+ 


3.3 2 . 

r y (m -3m -9m-5) 


16m 7 / 2 


— cosh -1 _ cosh - 
m+1 ' 


cosh 


— i ( 2m-ram b +m fc ) y+2mm fc c c 


HL{.(lIH-l)y 

. (2nL+ma^-Si ^ ) y+2mm^c 0 


m t (m+l)y 




y-HntCcriaty) (my4mmtc 0 +mty) 


B(y-Hntco) 2 
2i<m m^ 3 


_ 2mm t y(y4-m t c 0 ) (7m+5)-m t g y g C3m g +22m+15) 

2W 3 mt 3 


y 3 (3m 2 +8m+13) ^(l-m) 


2km 


T/2 


] 


xdx 
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The constant F 3 is evaluated as: 



which yields for mt / 1 


F 0 = 


- m 3 f 3 — (m 2 f g +2Pm 2 f— g g a) 3/2 p 3 


6xat 2 m 2 a 


6m" 


cosh - 


m 


(P+f) 

Pmt 


f (ag— 3m s f) ^mgfg+2pm2f-apg mff 

6m 2 a 2 2a 2 


f 3 (2m t 2 +m 2 ) 
+ 6a 5 / 2 



m 2 f-ap 
cos 

™*t 2 c Q 


mf \ 


and for m^ = 1 


- m3c o 3 -[ m2c o 2+2 P m2c cr-P 2 (l-^ 2 ) 1 3/2 , _pf_ m(P+c Q ) 

6m 2 ( 1-m 2 ) 6m 3 p 

c 0 [p(l-m 2 )— Sn^co] */n 2 c 0 2 +2pm 2 c 0 — P 2 (l-^n 2 ) _ • mc 0 3 
+ 6m 2 ( 1-m 2 ) 2 + 2(l-m 2 ) 2 


c 0 5 (24-m 2 ) 
• + 6(l-m 2 ) 5 / 2 


cos 


m 2 c 0 — P(l-m 2 ) 
mc 0 


cos 1 m 
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Also N4 Is evaluated as: 


r>3 

*V = / f4(y)^y 

J Ci’ 


which yields for mt f 1 


*4 = 


L24m 3 m^ 3 ^ 6h J l8m 3 m^. 3 J 


+ 2£\ + _L 


(mf g +3df—3 2 b) 3 ^ a 

-4b 


+ t 


r2fiMf + 5^?^ 


4b 


8b 3 / £ 


df— 2pb 

COB COB - 

mtf(m+l) 


m 1//2 df 2 1 

+ 4b } 


■ Af 2 ( 5d. 2 +4mb ) +SBbdf +128m a b 2 f 2 
384m a m t S b a " 


( 5Adf ‘ 

B \ 

[(mf 2 ) 3 / 2 " 

\ l44m 3 nL t - 3 b 

18 m 

4b 


(0 




cosh- 1 J!&2* - cosh- 1 iS - cosh" 1 

Pmj.(in+1) m+1 pm^m+l) 


3 2 5Pdf 5d g f g 2mf g 

12 b 2 8b 3 3b 2 


r £ 

\3b 




*/mf 2 +3 df — 3 2 b ) 




(ii + 5p + ^ + i f; 

df 3 ( 12mb— 5d 2 ) 




! df— 20b -1 d 

C0B j- — 1 cos — - - 

m^fCm+l) m^.(m+l) 


d 

mt(m+l) 
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hf‘ 3 (l2mg+5h. 2 ) 

l 6 g^ 



-1 hf— 2f3g 

cos" 

DL(.f(m4-l) 

2m 

_ 5^ , 2m 

3b 2 " 

Sg 3 Bg 2 / 


— L 


1 




Df 2 ( $h 2 +4mg) -fSHfgh+lgSf^g 2 | f/i 2pg-flA / vW^Bfhr-B 2 ® ) 

BSta^m^g 2 J L\ 4g J 


/f^tim 1 / 2 \ rm t 2 f 2 (m+l) 2 "I f fh-23g h ~| "1 

+ v~^r~/L "' "s? 575- J L cos ^i?(^iT “ C08 ^irJj 


Z' D 
[_ \ 24m 3 mt 3 . 


.)("*Hdb)] 


(mf g +3fh-p g g) 3 / a 

4 g 


e 4 (3gf+am+I3) ^(l-m) 1 . f (mf 2 ) 372 ^ Z 7 5Dfh YZ _JL_Y 1 

96m 7 / 2 J L 4g JL\ l44m 3 mt 3 g/ \ lSm^t 3 / J 


•where 


A = fika 2 +l4m 2 mt+10mmt+3m 2 mt 2 +22nmit 2 + 15 m t 2 
B = 2mffltc 0 (8m.+Tini]it + 5mt) 

D = 8m 2 -l4m 2 mt-10™t+3m 2 mt 2 +22nnnt 2 +15mt 2 
H = 2mm t c 0 (8mr-7ima^-^m t .) 
and for = 1 


F, 


= I— 

baa 3 


?c 0 ( 3 mr-l) 
12(mr-l) J 


. 7 1 r '[mco 2 +3c 0 (3mr-l)-f2P 2 (iiir-l) ] 3/2 1 

24m 3 J \ 6(nt~l) J 


r t 

5c 0 (3ni-l) 

7 1 

f(mc 0 2 ) 3 / 2 l 

+ f T 

c 0 2 ( 37m 2 -22m+5) 

L32m 3 

_ 12(m— l) 

24m 3 J 

6(mr-l) 

W 3 

64-(m-l) 2 
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V 

c 0 (3m-l) 

24m 3 

L 4(m— l) 


+ c 0 2 ~l J 4p (m— 1) + c 0 ( 3m-l) 

3m J l 8(m-l) 


*Jmc 0 2 +pc 0 ( 3m-l)+2P 2 (m— 1) 


j, c 0 2 (m+l) 2 
+ l6V2(l-m) 3 / 2 


, 4p(mr-l)+c 0 (3mr-l) 3m-l 

cos j- - cos 

c 0 (m+l) m+1 . 


m 1/2 c 0 2 (3m-l) 
8(l-m) 


■} 



"cosh- 1 _ cosh- 1 

3 (m+1) 


3~m _ cosh -1 P(3nj-l)+2mc 0 
m+1 p(m+l) 


+ (mc 0 ) 


2p 3 v^co 2 +pmc 0 (m+l) 


7mc 0 (m+l) 


35(l+m) 4 

r P 2 _ 5pco(3m-l) 
j_6m(m-l) 48m(m-l) 2 


+ • 


5c 0 g (3m-l) 2 
64m(m— l) 3 


Co 


6(mr-l) 2 J 


J 2p 2 m(m-l) +Pmc 0 ( 3m— 1) +m 2 c 0 2 


L 3£m 3 co 3 j_ c o 3 ( 63m 3 — 39m 2 +2Im-5 ) ^ nfl -i 4p(m-l)+c 0 (3m-l) 

35 (m+1) 4 128 V^(l^(mr-l) s L c D (m+l) 


— COB 


1 


3m-l 
m+1 _ 


5c 0 3 (3m— l) g 
64(mr-l) 3 



P 4 (3m 2 +8m+13) ^(l-m) ' f 2c 0 ^ (5-m) [2mco-3P(m+l) ] 

96m 7 / 2 ' _ 1 9m(m+l) 90m 2 (m+l) 2 
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(5-m)C8m 2 c 0 2 -123mc 0 (m+l)+15P 2 (m+l) 2 3 \/ r P „ , ,m 3 / 2 I 

1 i^o +p ° 0 ( "D 1 j 


(men 2 ) 3 / 2 ] f -g£o- - + C pfe sl. + ' 

J L9m(nH-l) 45m(m+l)2 105m(m+l)2 


■where 


T = 25m 2 +32ntfl5 
V = 10mc o (3m+l) 
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V 



Mach lines 


Figure /. — Geometric characteristics of the wing beam. 








Figure J. — The boundary conditions and 
solution for differential twist. 



pressure distribution for the conical flow 
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(a) Span toad distributions for the rigid and elastic wings 
at the same angle of attack. 



Percent semispan 


(b) Span toad distributions for the rigid and elastic wings 
at the same total lift coefficient . 

Figure 5. — A comparison of the span load distributions for the rigid 
and the elastic wings in accelerated flight. M, 1.414 } omax, 30,000 
pounds per square inch ; nW/S , 150 pounds per square foot } 
q, 2 II. 2 pounds per square foot. 
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(a) Mach number, 1.414. 



Figjre 7 — Variation with dynamic pressure of the ratio of the parameters C m 
for the elastic and rigid wings in accelerated flight. a 




Neutral point shift forward due to wing elasticity in Neutral point shift forward due to wing elasticity in 
accelerated flight, percent mean aerodynamic accelerated flight, percent mean aerodynamic 

chord, % a chord, % 7T 
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Dynamic pressure , q, hundreds of pounds per square foot 
(a) Mach number , 1.414 . 



(b) Incompressible flow. 

Figure 8. - Variation with dynamic pressure of the neutral point &ift 
for the elastic wing in accelerated flight. 



